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Abstract. Optical light curves of bright BL Lac objects, in particular those having the SED synchrotron peak at frequencies 
lower than 10'"* Hz, are characterized by long-term trends with time scales of a few decades. These variations are probabily not 
related to fast perturbations moving down the jet, which have a much faster evolution, but more likely to slow changes of other 
physical or geometrical parameters, for example the jet direction. In these cases we can expect changes of the inner jet structure 
detectable with VLBI imaging. Recent literature results are reviewed and discussed in this context. 
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1. Introduction 

BL Lacertae objects are AGNs characterized by a Spectral 
Energy Distribution (SED) having a non-thermal emission over 
the entire electromagnetic spectrum. The typical SED shows 
two broad peaks, one at lower frequencies generally explained 
by Synchrotron Radiation (SR) from relativistic electrons mov- 
ing along a jet closely aligned to the observer's line of sight 
(Blandford and Rees 1978), and one at higher frequencies orig- 
inated by Inverse Compton radiation (ICR) from the same elec- 
tron population. Padovani and Giommi (1995) classified BL 
Lac objects in two groups based on the maximum frequency 
vs of SR peak. The sources with vg in the range 10'^ - 10'^ 
Hz are called Low-Energy peaked BL Lacs (LBL), those with 
Vs ^ 10'^ - 10'^ Hz are named High-Energy peaked BL Lacs 
(HBL). 

LBL sources generally show large variations in their emis- 
sion over a wide range of time scales from minutes to years. 
Their fast variability, high apparent luminosity and a brigth- 
ness temperature exceeding the Compton limit (Kellermann 
and Pauliny-Toth 1969) imply relativistic boosting, defined by 
the Doppler beaming factor D - l/r(l - fi cosQ), where 
r = (1 - ffy^^^ is the bulk Lorentz factor of the electrons 
in the jet and 6 the angle to the observer's line of sight. A good 
estimate of D is relevant to understand the physics of BL Lacs. 
Furthermore, D can change during the source life time imply- 
ing large variations of the mean apparent flux (S oc D^^"^). 

Optical emission in LBL sources is thought to be the high 
frequency tail of the SR component. Therefore, one can ex- 
pect to observe a correlation with changes observed in the ra- 
dio band. When historic optical light curves are available, it is 
found that variations can also occur over time scales of several 
decades. Large flux changes with much shorter time scales are 
typically superposed to these long term trends. In particular, the 
secular variations can be reasonably associated with changes in 
the structure and/or direction of the inner jet. VLBI images are 
then very useful to obtain informations on the evolution of the 
Doppler factor D. 



In this contribution we present the possible correlation be- 
tween changes in the optical luminosity as observed for LBL 
sources and the variability observed in their VLBI radio struc- 
tures. Many of these sources are also important targets for 
the coming space missions for high energy astrophysics like 
AGILE and GLAS T. 



2. Optical historic light curves and monitoring 

Historic light curves of BL Lac objects can be obtained only 
from photographic plate archives. They are usually not well 
sampled over the maximum time interval of about a century. 
An example is the B light curve of OQ 530 shown in Fig. 
1 (Massaro et al. 2004a). Photographic magnitudes were ob- 
tained by Miller ( 1 978) from the Harvard archive and have been 
corrected to match the photometric B band (Nesci and Massaro 
1999). Note the long term trend with a mean decreasing rate of 
about 15-20% in ten years. 

For the well known source OJ 287 the analysis of the his- 
toric data indicated the existence of recurrent large flares with 
a period of 11.95 years (Sillanpaa et al. 1988) confirmed by 
the outburst observed in 1994. A photometric light curve of 
OJ 287 in the R, covering the past 30 years, is shown in Fig. 
2. Note, in addition to the main oubursts of 1982 and 1994, the 
decreasing trend up to Spring 1999 when the source luminosity 
reached its minimum. After that epoch the brightness of OJ 287 
started to increase at a rate higher than that of the decreasing 
part. These long term trends have been indicated in Fig. 2 by 
the dashed tick line drawn following the local minima. Note 
also that the highest flux level reached by OJ 287 in the out- 
burst of 1994 is about one third of that of 1982 maximum and 
seems roughly proportional to the mean flux. If the periodic re- 
currence of these oubursts will be confirmed in future, the next 
one, expected for 2006, could reach a flux comparable or even 
higher than that of 1982. 

Another well studied source is ON 231 (Massaro et al. 
2001). Its behaviour was characterised by an increase of the 
mean brightness after a minimum in the early seventies culmi- 



40 



E. Massaro and F. Mantovani: Optical and Radio emission from BLLac objects 




1900 1950 2000 



time (years) 

Fig. 1. The historic optical light curve of OQ 530 in the B band. Open circles are the data from the Harvard archive (Miller 1978) 
corrected to match nearly simultaneous data from Asiago archive (filled circles). The most recent points are photometric data 
from (Massaro et al. 2004a). 



nated in a very strong outburst in April-May 1998 (Massaro et 
al. 1999). After that event the mean source luminosity slowly 
declined, continuing to show changes greater than one magni- 
tude (Tosti et al. 2002). 

The amount of available data on variability of bright 
LBL objects has largely increased in the past ten years 
thanks to the use of small aperture telescopes equipped 
with CCD detectors. It is now simple to make these in- 
struments fully automatic obtaining a very high observa- 
tional efficency. Automatic optical monitoring provided very 
well sampled light curves for several sources (see for 
instance http://astro.fisica.unipg.it/PGblazar/tabella2000.htmi. 
They provide prompt information on the onset of active phases, 
important to trigger observations from space X- and y-ray ob- 
servatories. The best known case is that of BL Lac itself which, 
after a rather quiescent phase, in 1997 showed a very intense 
activity characterized by a very high luminosity and a large in- 
traday variability. 

In the radio band LBL objects are also variable but gen- 
erally not so violently variable as in the optical. Monitoring 
observations in the GHz range are important to search for 
correlations between optical and IR data, useful to investigate 
the radiative life times of relativistic electrons. The most 
continuous radio monitoring programs have been carried out 
for several years by the groups of the University of Michigan 
(http://www.astro.lsa.umich.edu/obs/radiotel/radiotel.html) 



and of the Metsahovi Radio Observatory 
( http : //kurp- www.hut.fi/quasar/ , Ciaramella et al. 2004). 

3. The Spectral Energy Distribution 

To better understand the link between the optical and the radio 
emission in LBL sources it is important to have a good de- 
scription of the broad band spectral properties of the SR. Non- 
thermal spectra are generally well fitted by single power laws 
over frequency intervals of the order of one decade. However, 
on wider intervals power law models fail to reproduce the ob- 
served mild curvature. Landau et al. (1986) showed that the 
spectra of some BL Lac objects from radio to the optical and 
UV ranges are represented by a parabola in lograrithmic coor- 
dinates 

5(v) =A(v/v„)-<"+''^"«(''/''''» (1) 

better than an exponential cut-off. Here the SED curvature is 
described only by the parameter b. This spectral law has been 
used successfully to fit the X-ray SEDs for the HBL sources 
Mkn 421 and Mkn 501 (Massaro et al. 2004b,c) in a variety 
of luminosity states. This implies that the electron spectra also 
show a similar distribution, which can be obtained from statis- 
tical acceleration mechanisms when the probability for a parti- 
cle to remain inside the acceleration region decreases with the 
energy of the particle itself. Alternatively, when the particle en- 
ergy remains below a critical value the acceleration probability 
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Fig. 2. The photometric R light curve of OJ 287 in the R band (not corrected for the interstellar reddening). The two outbursts 
of 1982 (around JD= 5,400) and 1994 (around JD=9,600) are evident. Note the mean decresing trend (dashed line) down to the 
minimum observed in March 1999 followed by a faster increase continuing up to now. (Thanks to N. Marchili and G. Tosti). 



is constant and therefore the resulting spectrum is given by a 
power law and changes into a log-parabola when the accler- 
ation probabilty decreases. The resulting SR spectrum can be 
then approximated by 

5(v) = Av"^°+''^°«''^''^''''^ . (2) 

Some spectra of the SR observed from LBL objects are 
plotted in the upper panel of Fig. 3. Equation (2) has been used 
for the two spectra of OQ 530, while Equation (1) for the re- 
maining sources. 

Another interesting possibility is that presented in Fig. 3 
(lower panel) in which there are two SR components peaking 
at different frequencies. The two spectra were computed using 
the same values of b as the electron populations were acceler- 
ated by similar processes. They could be originated in differ- 
ent regions along the jet and evolve with different time scales. 
Simultaneous broad band observations may not provide enough 
information to distinguish the two components, while variabil- 
ity studies can give the most useful information. 

4. VLBI imaging 

VLBl imaging of BL Lac objects is important because we can 
obtain direct information on the evolution of the inner jet struc- 
ture. Furthermore, polarisation measurements allow us to know 
the orientation of the magnetic field (see for instance Gabuzda 
et al. 1999, 2004). At the typical distances of bright BL Lacs 



objects, say z < 0.3, the VLBI angular resolution at 5 GHz 
corresponds to distances of few parsecs. As a consequence, 
only structure changes involving time scales of years can be de- 
tected. Images at higher frequencies (say 43 GHz) are very use- 
ful to investigate the variations on shorter time scales. However, 
observations at lower frequencies are needed to define the spec- 
tral behaviour. Savolainen et al. (2002) applied a model for 
the time evolution of radio flares of a blazar sample, including 
some LBL objects, observed at 22 and 43 GHz. Their analy- 
sis was successful in describing events with typical time scales 
of one-two years. They showed that in more than 50% of the 
sources, a new VLBI component appeared in the jet after a to- 
tal flux density flare. Usually these components brighten and 
decay within ~0.15 mas of the radio core. This result supports 
the interpretation that flares are originated by a sudden energy 
injection from the core. Subsequent shocks propagate in the jet 
and are associated with superluminal components. For faster 
variations, like those frequently observed in the optical for LBL 
sources which evolve on time scale of a month or even shorter, 
new components emerging from the core cannot be easily seen. 
They are likely exhausted before they can reach an angular sep- 
aration large enough to be resolved out by VLBI. 

Long term optical trends are likely originated by different 
and more regular processes. An interesting possibility is that of 
a precessing jet, which implies a change of 6 and consequently 
of the Doppler factor D. One can expect to observe a relation 
between the superluminal velocity of the new born components 
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Fig. 3. Upper panel: Spectral Energy Distribution of some LBL 
objects, (adapted from Perri et al. 2003 and Massaro et al. 
2004a). Lower panel: a possible model of a two component 
SED. 

and the observed total core luminosity. This relation is shown 
in Fig. 4 where the Doppler factor is plotted as a function of 
Papp (Blandford and Konigl 1979). For smaller than about 8°, 
an increase of 6 would imply a higher superluminal velocity 
corresponding to a fading of the luminosity and viceversa. Of 
course, it is not easy to separate such an effect from those due to 
short time changes. Dense radio-optical and VLBI monitorings 
(typically twice a year) are here required. 

Evidence for a precessing jet nozzle in BL Lacertae has 
been recently found by Stirling et al. (2003) using a series of 
VLBA high resolution and polarisation images at 43 GHz taken 
from March 1998 to April 2001. They estimated a period of 
~2 years. There is a very large amount of optical data on this 
source, in particular after the large outburst of 1997, however, 
a complete analysis has not been carried out. Fan et al. (1998) 
analysed the historic light curve of BL Lacertae and found sev- 
eral recurrence times (one of them of about 2 years). Further 
investigations to confirm this result, in particular on the phase 
of the periodic changes, are necessary. Furthermore, this anal- 
ysis would be useful to understand if the spectrum of the SR 
observed in the optical extends to the radio band or if it is gen- 
erated by a different component. 

EVN images at 5GHz of OQ 530, recently presented by 
Massaro et al. (2004a), show a bright radio core and a rather 
weak jet. Optical and radio light curves are charaterized by a 




Fig. 4. The Doppler factor as a function of the apparent velocity 
(in units of the speed of light) for = 2°, 5 MOM 5° (from top 
to bottom). The other lines show the relation between D and 
Papp when is changing and F is taken equal to 5, 8 and 10. 

low luminosity phase at the beginning of 1999 and they start to 
brighten later on. Simultaneous optical and radio observations 
(February 1999 and June 2001) indicate a similar increase of 
the core flux in the R band and at 5 GHz (Fig. 3 in Massaro et 
al. 2004a) consistent with a unique SR emission component. 

An interesting point to rise discussing the structure of BL 
Lac objects (and more in general for blazars) is the presence of 
stationary components detected near the core. In the extensive 
VLBI monitoring at 22 and 43 GHz on a sample of 42 y-ray 
loud blazars, including 11 LBL objects Jorstad et al. (2001) 
and Marscher et al. (2002) found stationary components within 
a distance of 2 mas from the core in 27 sources. They proposed 
that those components could be associated with a standing hy- 
drodynamical compression. In such a case some components 
may appear blended in a single feature without a significant 
proper motion making it diflicult to detect superluminal mo- 
tion in the inner part of the jet. 

Finally it is worth to mention the case of ON 23 1 (Massaro 
et al. 2001; Mantovani & Massaro, this conference) were a 
component not previously detected was observed on the op- 
posite side of the jet with respect to the radio core. Again, the 
change of the direction of the jet, occurring after the Spring 
1998 optical flare, was considered as a possible explanation for 
such unusual structure. 

5. Conclusions 

The next ten years will be very interesting for the astrophysics 
of BL Lac objects. These sources and Flat Spectrum Radio 
Quasars, grouped together in the Blazar class, are charaterized 
by a strong and variable non-thermal emission. For this rea- 
son they are the most numerous extragalactic y-ray sources and 
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are also well detected in the IR range. The majority of the 208 
sources detected by WMAP (Bennett et al. 2003) are indeed 
blazars (Giommi & Colafrancesco 2004). 

Two high sensitivity space observatories, well suited to ob- 
serve BL Lac objects, will be operative in the next three years: 
GLAST for y rays above -30 MeV and Planck (30 GHz - 
857 GHz). It is expected that at least a few thousands sources 
will be detected. It is important at this stage to select a sample 
of target sources planning a monitoring programme both in the 
radio and optical bands. In particular, it will be very useful to 
plan for a monitoring programme with high resolution radio in- 
terferometers to obtain detailed information on milli-arcsecond 
scale structural changes and brightness variations of the broad- 
band synchrotron emission. 

Optical variability plus VLBI imaging can also be useful to 
verify the unified scheme model for BL Lac objects. According 
to that interpretation, the relativistic effects are relevant when 
the angles between the jet direction and the observer's line 
of sight are small. When this angle increases the source be- 
comes fainter, its variability time scale becomes longer and, 
at large angles the source should appear as a FRI radiogalaxy 
(Ghisellini et al. 1993). Some BL Lac objects are rather faint 
in the optical while their radio flux is relatively high. The syn- 
chrotron peak frequency of these sources is rather low, < 
lO'^ Hz. We will refer to them as Very Low energy peaked BL 
Lac objects, (shortly VLBL). In the unification scenario VLBL 
objects could have Doppler factors smaller than that of typi- 
cal LBL sources. Their behaviour should be intermediate be- 
tween these sources and FRI radio galaxies. Unfortunately, too 
few data and images are available for this 'class' of objects. 
Therefore no firm conclusions can be at present derived. It wiU 
be then useful to plan coordinated observational programs for 
a sample of VLBL objects to investigate their main properties. 
In particular, simultaneous broad band observations are needed 
to have information on the SED and VLBI observations the de- 
tect superluminal motions of the jet components to estimate the 
Doppler factor. 
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